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PREFACE 

Part  of  the  RAND  research  program  consists  of  basic 
supporting  studies  In  mathematics.  The  mathematical 
research  presented  here  concerns  techniques  for  the 
solution  of  differential  equations  with  time-lags. 

This  method  is  of  Importance  In  connection  with  the  study 
of  more  realistic  models  of  chemotherapy,  of  the  type 
being  studied  under  GM-09608. 


SUMMARY 


In  this  paper  we  briefly  indicate  how  a  technique 
for  the  reduction  of  the  solution  of  differential- 
difference  equations  with  one  time-lag  to  the  solution 
of  systems  of  ordinary  differential  equations  can  be 
extended  to  the  more  complex  situation  involving 
different  time-lags. 
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A  NOTE  ON  THE  COMPUTATIONAL  SOLUTION  OF  A  SYSTEM  OP 
DIFFERENTIAL  EQUATIONS  WITH  VARYING  TIME-LAGS 

1 .  INTRODUCTION 

In  this  paper  we  wish  briefly  to  Indicate  how  a 
technique  for  the  reduction  of  the  solution  of 
differential— difference  equations  vrlth  one  time-lag  to 
the  solution  of  systems  of  ordinary  differential 
equations  can  be  extended  to  the  more  complex  situation 
involving  different  time-lags.  This  method  Is  of  some 
Importance  In  connection  with  the  study  of  more  realistic 
models  of  chemotherapy,  of  the  type  studied  in  [  5]  ,  [  6]  . 

2.  DESCRIPTION  OF  METHOD 

Consider  a  system  of  tvro  differential— difference 
equations  with  different,  but  commensurable,  time-lags, 
as  and  bs,  where  s  >  0  and  a  and  b  are  positive 
Integers: 

(2.1)  y'(t)  =  f(t,y(t),y(t  -  as),z(t),z(t  -  bs)), 

z'(t)  «=  g(t,y(t),y(t  -  as),z(t),z(t  -  bs)), 

with  the  Initial  conditions  y(t)  =  z(t)  =0,  t  ^  0. 

Let  us  now  present  an  algorithm  v^lch  obtains  the 
nxmierlcal  solution  of  (2.1)  by  means  of  the  solution  of 
sequences  of  ordinary  differential  equations  of 
Increasing  order.  The  case  of  one  time-lag  was  treated 
In  [1] ,  [2]  . 
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Int reduce  the  functions 

(2.2) 

Zj^(t)  =  z(t  +  ns), 

with  n*=  —  b,  —  b  +1,  —  b  +  2,...,0,  ...,N,  where  N 
is  determined  by  the  range  of  integration,  and 
0  ^  t  s. 

To  begin  with,  we  solve  ninnerically  the  following 
system  of  ordinary  differential  equations: 

(2-3)  yJjj  =  0,  =  0, 

y-b+l  “  °* 

« 

yli  = 

^6  “  s(t,yQ,y_^,ZQ,z_^ 

in  the  range  0  ^  t  £  s,  with  the  initial  values 
yj^(O)  =  ZjL(O)  for  all  1. 

Note  that 

(2. 4)  yi(0)  =  yQ(5),  z^(0)  =  Zq(s). 

In  the  next  roxmd,  we  adjoin  to  the  system  in  (2.3) 


the  eqviatlons 
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(2.5)  “  f(t  +  s>yi*y_a+ii* 

yi(o)  -  yQ(s), 

*i  -  g(t  + 

Zj^Co)  -  Zq(s). 

Integrating  this  larger  system  over  the  Interval  [Ojs], 
we  obtain  two  more  Initial  values 

(2.6)  y2(0)  =  Yj^Ca),  z^io)  -  z^(s). 

Generally,  In  the  (1  +  1)— st  round  of  Iteration,  we 
add  to  the  1— th  system  the  equations 

(2.7)  yj^  «  f(t  +  yi(o)  «  yi_i(s), 

z|  -  f(t  +  i8,y^,y__^^j^,z^,z_^j^^),  Zj^(o)  »  z^_j^(s). 
3.  DISCUSSION 

In  many  cases  the  time-lags,  and  Xg,  will  not 

be  commensurable.  If  we  do  not  wish  to  use  the  usual 
storage  technique  for  the  numerical  Integration  of  (2.1), 
we  can  approximate  to  X^  and  Xg  by  rational  fractions 
Pi/q*  P2/<1  such  a  way  that 

Pi  2  p«  2 

(3.1)  (X^  -  -^)  +  (Xg  -  -|) 

Is  as  small  as  possible,  with  q  ^  N.  If  q  Is  too 
lao’ge,  the  systems  of  equations  described  In  Sec.  2 
become  unwieldy. 
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In  general,  we  can  cut  down  on  the  dimension  of 
these  differential  equations  by  applying  the  preceding 
technique  over  an  interval  [  and  then  using  the 

technique  of  differential  approximation  [ 3i ,  [  to 
replace  storage  of  the  "initial  values"  over  [t^,T^  +  si. 
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